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We consider topological invariant describing the vacuum states of superfluid 3 He-B, which belongs to the 
special class of time-reversal invariant topological insulators and superfluids. Discrete symmetries important 
for classification of the topologically distinct vacuum states are discussed. One of them leads to the additional 
subclasses of 3 He-B states and is responsible for the finite density of states of Majorana fermions living at some 
interfaces between the bulk states. Integer valued topological invariant is expressed in terms of the Green's 
function, which allows us to consider systems with interaction. 



PACS: 

1. INTRODUCTION 

General classification schemes based on topology 
[2 [3j H OH El [7] suggest existence of the topolog- 
ical insulators and fully gapped topological superflu- 
ids/superconductors which have the gapless edge states 
on the boundary or at the interface. Superfluid 3 He- 
B belongs to the special class of three-dimensional 
topological superfluids with time-reversal symmetry. 
Topological invariant which describes the ground states 
(vacua) of 3 He-B has been discussed in Refs. [1] and [§]. 
Here we present the explicit expression for the relevant 
topological invariant, and discuss topological quantum 
phase transitions occurring between the vacuum states 
and fcrmion zero modes living at the interfaces between 
the bulk states. 



2. TOPOLOGICAL INVARIANT PROTECTED 
BY SYMMETRY 

Usually in the topological classification of ground 
states people use the Hamiltonian of free particles or the 
corresponding effective Hamiltonian such as Dirac and 
Bogoliubov-de Gennes Hamiltonians [TJ [31 [3J 13] . How- 
ever, in this classification the natural problem arises, 
what is the effect of interaction between particles. More- 
over, the original first- principle many-body Hamiltonian 
of, say, liquid 3 He 
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has no information on the topological structure of 
the ground state of the system - superfluid 3 Hc-B. 
The accurate procedure to reduce such a strongly in- 
teracting many-body system to the effective coarse- 
grained Hamiltonian is absent. However, the micro- 
scopic Hamiltonian allows us at list in principle to calcu- 
late the Green's function G(oj, p) - the quantity which 
determines the main properties of the translational in- 
variant or periodic ground states of the system. That 
is why the object for the topological classification must 
be the Green's function rather than Hamiltonian. Then 
it is applicable even in cases when one cannot introduce 
the effective low energy Hamiltonian, for example when 
Green's function does not have poles, see [HE]. 

Green's function topology has been used in particu- 
lar for classification of topologically protected nodes in 
the quasiparticle energy spectrum of systems of differ- 
ent dimensions; for the classification of the topological 
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ground states in the fully gapped 2 + 1 systems, which 
experience intrinsic quantum Hall and spin-Hall effects 
[ini HU [HI [13j, including multi-band topological 
insulators [141 : in relativistic quantum field theory of 
2 + 1 massive Dirac fermions [HI [HI [18] ; etc. 

The integer-valued topological invariants were ex- 
pressed via the Green's function, which was considered 
at imaginary frequency to avoid the zeroes, poles or 
other possible singularities on the mass shell. In the 
fully gapped systems one may consider the Green's func- 
tion not only at imaginary frequency, but also at real 
frequency if it is below the gap. For the classification 
of the 3 He-B states we shall use the Green's function at 
zero frequency: G(p) = G(cu — 0,p). The typical ex- 
ample of the integer valued topological invariant is the 
following 3-form: 



N 



24tt 2 



tr 



d 3 P Qd Pi g- l gd Pj g- l Gd Pk g 



(2) 



where integration is over the whole momentum space 
for translational invariant systems, or over the Brillouin 
zone in crystals. For 3 He-B obeying time-reversal sym- 
metry this invariant is identically zero, N = 0. How- 
ever, the discrete symmetries of 3 He-B give rise to the 
other invariants. Examples of additional integer valued 
topological invariants which appear due to symmetry 
in different condensed matter systems and in quantum 
field theory can be found in Refs. [H [H QH G3 1] • 

Due to symmetry, one or several Pauli matrices 7 
of spin, pseudospin or Bogoliubov-Nambu spin may ei- 
ther commute or anti-commute with the Green's func- 
tion matrix: 

yQ(p) = -S(ph , (3) 



or 



7<?(p) = Q(ph- 



(4) 



This leads to the following integer valued topological 
invariants: 



N lilLtr 
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■Jd 3 P gd^g-'gdp^gdp.g- 1 
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with g(p) and 7 obeying either §5§ or (H]). 

The topological classes of the 3 He-B states can be 
represented by the following Green's function, which 
plays the role of effective Hamiltonian: 



L (p) = M(p)r 3 + n (a x c x p x + (JyCyPy + a z c z p z ) , (6) 

^(P) = ^r-M,(7) 



2m* 

where n are Pauli matrices of Bogolyubov-Nambu spin. 
The overall 'conformaP factor, which may depend on p, 



is omitted since it docs not influence the topological in- 
variant. In the isotropic 3 He-B all 'speeds of light' are 
equal, \c x \ — \c y \ = \c z \ = c. The topological invariant 
relevant for 3 He-B is iV 7 in with 7 = T2: 



N 7 = 



24tt 2 



tr 



t 2 J d 3 P gd p ,g- x gd Vj g- l gd Vh g- 



(8) 

The Ti matrix plays the role of the 7-matrix, which 
anti-commutes with the Green's function in (|3|). The T2 
symmetry is combination of time reversal and particle- 
hole symmetries in 3 He-B. 



3. PHASE DIAGRAM IN MASS PLANE 

Fig. Q] shows the phase diagram of topological 
states of 3 He-B in the plane [p,, 1/m*). On the line 
1 /to* = one obtains Dirac fermions with mass param- 
eter M = —fj,: 

g^ 1 (p) = Mt 3 + n (a x c x p x + (J y CyPy + a z c z p z ) . (9) 

Topological invariant ijHJ) for Dirac fermions is 



7V 7 = -sign(Mc x CyC ZJ 



sign(fic x CyC z ) . 



(10) 



In relativistic quantum field theory ti matrix in (jSj) 
is the operator of CT symmetry [50]. Hamiltonian for 
Dirac fermions is non-compact, with different asymp- 
totes at p — > 00 for different directions of momentum 
p. As a result, the topological charge of the Dirac 
fermions is intermediate between charges of the com- 
pact states of 3 He-B below and above the horizontal 
axis (see Refs. [21], p] and [5] on the marginal be- 
havior of fermions with relativistic spectrum; also note 
that the topological invariant iV 7 in © has values twice 
larger than the invariants introduced in Refs. [3] and 
[5]). When the line 1/m* = is crossed, the asymptotic 
behavior of the 3 He-B Green's function changes from 
£ _1 (P) -» +T 3 p 2 /2m* to g-^p) -» -T 3 p 2 /2m*. 

The real superfluid 3 He-B lives on the weak coupling 
side of the phase diagram: at fi > 0, to* >0, /i> m*c 2 . 
However, in the ultracold Fermi gases with triplet pair- 
ing the strong coupling limit is possible near the Fes- 
hbach resonance [22- When /j, crosses zero the topo- 
logical quantum phase transition occurs, at which the 
topological charge iV 7 changes from A^ 7 = 2 to iV 7 = 0. 
Previously the topological quantum phase transition oc- 
curring at /J, = has been considered for the chiral 
p-w&ve states of the 3 He-A type: the transition from 
the gapless state at /i > to the fully gapped state at 
fi < [TH1[52]. In the 3 He-B case, both states are fully 
gapped, while the intermediate state at [i = is gapless: 
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Phase diagram of topological states of 3 He-B in equa- 
tion §6$ in the plane (/x, 1/m*) for the speeds of light 
c x > 0, c y > and c 2 > 0. States on the line 
1/m* = correspond to the Dirac vacua, which Hamil- 
tonian is non-compact. Topological charge of the Dirac 
fermions is intermediate between charges of compact 
3 He-B states. The line 1/m* = separates the states 
with different asymptotic behavior of the Hamiltonian 
at infinity: <5 _1 (p) — > ±T3p 2 /2m*. The line \i — 
marks topological quantum phase transition, which oc- 
curs between the weak coupling 3 He-B (with p > 0, 
m* > and topological charge N 1 — 2) and the strong 
coupling 3 He-B (with fi < 0, m* > and iV 7 = 0). This 
transition is topologically equivalent to quantum phase 
transition between Dirac vacua with opposite mass pa- 
rameter M = ±|/i |, which occurs when p crosses zero 
along the line 1/m* = 0. The interface which separates 
two states contains single Majorana fermion in case of 
3 He-B, and single chiral fermion in case of relativistic 
quantum fields. Difference in the nature of the fermions 
is that in Bogoliubov-de Gennes system the components 
of spinor are related by complex conjugation. This re- 
duces the number of degrees of freedom compared to 
Dirac case. 



This analog of the index theorem [23) implies that such 
interface contains single Majorana fermion. 

The same quantum phase transition occurs when fj, 
crosses zero along the line 1/m* = 0, i.e. the line of 
the relativistic Dirac fermions. At this transition the 
mass M = — [i of Dirac fermions changes sign. The 
rules (fTTj) and (fT2|) are also applicable for the Dirac 
fermions. However, in relativistic quantum field the- 
\1 ory the domain wall separating vacua with opposite M 
contains chiral fermion rather than the Majorana one. 
Difference in the nature of the fermions comes from the 
observation that in the Bogoliubov-de Gennes system 
the components of a spinor are related by complex con- 
jugation. This reduces the number of fermionic degrees 
of freedom compared to the relativistic quantum field 
theory and makes them the Majorana fermions. 

The spectrum of fermion zero modes at this interface 
can be easily found in the limit <C m*c 2 , when the 
p 2 term in the Hamiltonian can be neglected and the 
Hamiltonian approaches the Dirac one. In this limit 
Majorana fermions in 3 He-B have relativistic spectrum: 



H zm = cz ■ (o- x p) 



(13) 



The same spectrum is obtained for the chiral fermions 
living on the interface between the states of the 3D topo- 
logical insulator with normal and inverted arrangements 
of two bands, see Ref. [24] . 

The invariant ([5]) is applicable to the general case of 
systems with Green's function of the type 



Q'Hp) = e(p)r 3 + A(p)n 



(14) 



where e and A are Hermitian matrices with spin and 
band indices. This includes in particular the Hamilto- 
nians discussed in Refs. [3] and [25]. The phase diagram 
of topological 3D states in [3] is similar to that in Fig. Q] 
(note again that the topological invariant A 7 has values 
twice larger than invariants introduced in Refs. [3] and 
0). 



it has two point nodes at p = with opposite chiralities. 
There is a general relation between topological invari- 
ants of the two states and the number of point nodes 
in the intermediate gapless state [5], which for a given 
case reads: 

Appoint nodes = |A 7 (ju > 0) - A 7 (^ < 0)| = 2. (11) 

This relation also determines the number of 2+1 fermion 
zero modes living at the interface between the two 
states: 

N FZM = ^\N 1 (n 1 >0)-N 7 (u.< 0)| = 1. (12) 



4. PHASE DIAGRAM IN THE SPEED OF 
LIGHT PLANE 

Fig. [2] shows phase diagram of topological states of 
3 Hc-B in the plane (c x ,c y ) at fixed c 2 > 0, p > and 
m* > 0. When one of the components of speed of light, 
say, c x crosses zero, two pairs of point nodes appear in 
the intermediate gapless state, at points p = ±(0, 0,Pf), 
where p 2 F = 2pm* is Fermi momentum. This is the con- 
sequence of different topological invariant A 7 of the two 
states: 

Apoint nodes 

= |AT 7 ( C:c >0)-A 7 ( C;c <0)|=4. (15) 
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Phase diagram of 3 He-B states at fixed c z > 0, /i > 
and m* > 0. The interface between the states with dif- 
ferent winding number iV 7 contains Majorana fermions. 
However, in the presence of discrete symmetry Zi - 
combined rotation by it about z axis in spin and or- 
bital spaces - the 3 He-B states with the same topologi- 
cal invariant JV 7 may be connected only via the gapless 
state. As a result, the interface between two states with 
the same iV 7 also contains Majorana fermions. These 
fermions have finite density of states at zero energy [8]. 
This is the origin of the finite density of states of the 
Majorana edge states on the diffusive wall |27| . 



The intermediate state is the so called planar state [25] , 
and the corresponding gapless fermions in this state are 
characterized by the topological invariant protected by 
discrete symmetry in (U) 7 = T^a x , which commutes 
with Hamiltonian: jG(p) — S(p)7 (see Eq.(14.4) of 
Ref. 0): 



3 24tt 2 



r-.-.n, l dS 1 Gd p ^ 1 Gd Pv G^ 1 Gd Px G- 1 

(16)' 

where integral is around a nodal point p^ = (0, 0,pfj 0) 
01 Pn — (0,0, — Pf,0) in the 4D momentum-frequency 
space p^ = (p, w). Number of Majorana fermion zero 
modes living at the interface between the two states is 
correspondingly twice smaller: -/Vfzm = 2. The same 
number of 2 + 1 fermions live at the boundary of 3 Hc-B 
with specular reflection. This is because there is an 
exact mapping between the edge states on the 3 He- 
B boundary with perfect reflection and fermion zero 
modes living at the interface, see Ref. [27] , 



There are subclasses of 3 He-B states, which exist due 
to discrete symmetry. When two components of speed 
of light, say, c x and c y cross zero, the topological charge 
N 1 does not change. The final state can be continuously 
connected with original state by spin rotation by angle it 
about axis z. However, there is additional discrete sym- 
metry - the combined rotation by 7r about z axis in spin 
and orbital spaces. If this Z% symmetry is maintained, 
the two states with the same topological invariant 7V 7 
cannot be connected by adiabatic deformations. The 
intermediate gapless states necessarily contains the line 
of nodes. 

As a result the interface between states with the 
same iV 7 , at which two components of speed of light 
in (|10jl change sign, also contains zero modes. Spec- 
trum of these fermions H zm cx ca x k y , and they have 
finite density of states at zero energy [H [27] . Because of 
the mapping between the fermion zero modes living at 
this interface and the edge state of 3 He-B on diffusive 
boundary, the Majorana states on the diffusive wall also 
have finite density of states [27] in agreement with the 
detailed quasiclassical theory [28] . 



5. DISCUSSION 

There are many experimental evidences for the An- 
dreev surface states on the wall of 3 He-B, see e.g. 
[2"5l I5TJI |3"TI 152"] , and some other experiments are pro- 
posed [33]. At low temperature, where thermal quas- 
particles in bulk 3 He-B are frozen out exponentially, 
the low-energy surface states - Majorana fermions - will 
give the main contribution to thermodynamics and dis- 
sipation, with the power-law dependence of the physical 
quantities. The extra discrete symmetry of 3 He-B essen- 
tially enhances the effect of Andreev-Majorana states by 
providing the finite density of states of these modes on 
rough walls [34] . Systems with highly developed sur- 
faces, such as 3 He-B confined in porous materials, or 
systems of parallel plates or tubes, would be more ad- 
vantageous for the direct identification of the Majorana 
fermions in superfluid 3 He-B. 

It is a pleasure to thank A. Kitaev and Yu. 
Makhlin for valuable discussions. This work is sup- 
ported in part by the Academy of Finland, Centers 
of Excellence Program 2006-2011, the Russian Foun- 
dation for Basic Research (grant 06-02-16002-a), and 
the Khalatnikov-Starobinsky leading scientific school 
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